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Q ■ Abstract 

We develop quenched chiral perturbation theory for vector mesons made of 
light quarks, in the limit where the vector meson masses are much larger 
' than the pion mass. We use this theory to extract the leading nonanalytic 

■ dependence of the vector meson masses on the masses of the light quarks. By 

o 



ix. 



X 



comparing with analogous quantities computed in ordinary chiral perturba- 
tion theory, we estimate the size of quenching effects, observing that in general 
they can be quite large. This estimate is relevant to lattice simulations, where 
the p mass is often used to set the lattice spacing. 
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I. INTRODUCTION 



Lattice QCD has made great progress in recent years, to the point that the computation 
of the masses of the hght hadrons with the accuracy of a few percent is within sight [|l|. 
Of course, such a precision requires that one understand well the various approximations, 
extrapolations and ansatze which underlie the calculation, in order to estimate, and then 
minimize, the error which they induce. Without a doubt, the most significant of these, and 
by far the least well understood, are the corrections due to quenching, particularly as they 
affect the extrapolation of results to the physical values of the light quark masses. 

Quenching errors divide roughly into short and long distance parts. At short distance, 
their primary effect is to change the running of the coupling constant as- This manifests 
itself, for instance, in the wave functions of quarkonia, where it can be studied in detail 0. 
The long distance effects are by their nature more difficult to quantify. One approach 
is to use quenched chiral perturbation theory (QChPT), that is, the low energy effective 
theory of quenched QCD. This idea was proposed by Sharpe with the formalism 

further developed by Bernard and Golterman The advantage of using a version of 

chiral perturbation theory to study this problem is that its predictions follow from the basic 
properties of the underlying theory. Moreover, ordinary chiral perturbation theory is known 
to describe the low energy dynamics of unquenched QCD very well. 

More recently, QChPT has been extended to describe the interactions of pions with 
baryons and heavy mesons JTU|pT| in quenched QCD. In this paper, we will extend it 
further to describe the interactions of the light vector mesons, p, K* uj and (f). Jenkins, 
Manohar and Wise have formulated ordinary chiral perturbation theory for the vector 
mesons, treating the vector mesons as heavy particles, in a manner similar to chiral pertur- 
bation theory for baryons UTSf. We will adapt their formalism to quenched QCD. Although 



we will adhere to most of their conventions, we will use a slightly modified Lagrangian which 
generalizes more easily to the quenched case and to an arbitrary number Nj of light quarks. 

What can we expect to learn from such an investigation? The difficulty with low energy 
effective theories is that typically they contain a large number of undetermined coupling 
constants, all of which must be fixed from experiment. In fact, only for that sector of chiral 
perturbation theory which describes the self-interactions of the pions are there enough data 
to carry out such a phenomenological program beyond leading order. For all other particles 
and interactions in the chiral Lagrangian, the data are scanty at best. The situation in the 
quenched theory is even worse, since the only "data" available are extracted from lattice 
simulations. Hence we are left in a situation where firm and accurate predictions are largely 
impossible to obtain. 

Instead, we take a more modest approach. First, we note that there are certain loop 
corrections which are predicted unambiguously by the lowest order Lagrangian, namely 
those that have a nonanalytic dependence on the light quark masses. Such contributions 
to physical observables cannot be compensated by higher order terms in the chiral expan- 
sion. Furthermore, this nonanalytic dependence is particularly sensitive to the long distance 
behavior of the theory, which is what we wish to investigate. Second, we use the power 
counting rules which are implicit in chiral perturbation theory, by which all nonperturbative 
coupling constants are taken to be of order one (given certain other conventions). With 
this assumption, we can ask whether the nonanalytic dependence of a physical quantity on 
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the light quark masses is small or large. If it is small, we enhance our confidence that the 
low energy theory is well behaved. However, if it is large, then the applicability of the low 
energy theory to that quantity is called into question. 

The purpose of this investigation, then, is to explore whether QChPT is under reasonable 
control, by studying those quantities which the theory, at some level, predicts unambiguously. 
We will find, in fact, that QChPT is quite badly behaved for the masses of the vector mesons, 
for apparently reasonable values of the nonperturbative coupling constants. Whether this 
conclusion ought to be extended to quenched QCD itself is not completely clear, but at 
the very least, our results should sound a serious note of caution about quenched lattice 
simulations of the vector meson masses. 



II. QUENCHED CHIRAL PERTURBATION THEORY 

The extension of chiral perturbation theory to describe quenched QCD is by now stan- 
dard. Following Morel for each quark g*, a corresponding bosonic "ghost" quark is 
introduced with identical mass, so that the fermion determinant is canceled by the ghost 
determinant. These ghost quarks will then form mesons with the true quarks and with them- 
selves. As a result, the symmetry group of the theory is enlarged from SU{Nf) l x SU {Nf)ji 
to the semi-direct product {SU{Nf\Nf)L x SU{Nf\Nf)R) ® f/(l). 

In this larger theory, the matrix of Goldstone bosons is promoted to a supermatrix, 

n=(; i), (2.1) 

where the quark/ghost content of the fields is vr ~ gg, ^ QQ, X ^ QQ and tt ~ gg. Each of 
these is actually an Nf x Nf matrix; for example, for Nf = 3, vr is the ordinary pseudoscalar 
nonet 



TT 



^ V'h. (2.2) 



Note that x ^iid x^ are fermionic fields, while vr and vr are bosonic. 
The Lagrangian of QChPT is given by 

^Qx = Y (Str[5^Sa'^St] + ABo StT[M+]) 

+ ^ (Ao Str[a^n] Stiid^U] - Ml Str[n] Str[n]) . (2.3) 
Here S = ^ = e^^/^ (the normalization is such that /tt ~ 130 MeV) while 

M^C' I), (2.4) 



V M, 

M = diag(mi, . . . , niNf) , (2.5) 

and Ai± = ^{^^Ai^^ ± 'C-^O- The "supertrace" Str is defined with a minus sign for the 
ghost-antighost fields. The chief difference between the quenched and unquenched theories 
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1 /2 ~ 

is the presence of the terms involving Str[n] = A^^ {rj' — fj'). In the unquenched theory they 
can be neglected because they describe the dynamics of the rj' meson, which decouples from 
the theory. But quenching prevents the r]' from becoming heavy and decoupling, so these 
terms must be retained in QChPT. We have normalized Aq and Mq so that they have no 
implicit dependence on Nf. 

The propagators that are derived from this Lagrangian are the ordinary ones, except for 
the flavor-neutral mesons, where the non-decoupling of Ho leads to a curious double-pole 
structure. It is convenient to adopt a basis for the these mesons corresponding to qiq^ and 
qiQ^. Then the propagator in the flavor- neutral sector takes the form 



m 



where m\ = 2Bomi, and = 1 if i corresponds to a quark and = — 1 if i corresponds 
to a ghost. The second term in the propagator is treated as a new vertex, the so-called 
"hairpin," with the rule that it can be inserted only once on a given line. Note that this 
term can induce mixing between quark- ant iquark and ghost-antighost pairs. 

The treatment for vector mesons is similar. They are combined into the supermatrix 

with Hfj^ the usual matrix of vector mesons, which for Nf = 3 is 

^ \ 

1 









K*+ 

























^ Vs • (2.8) 

It is convenient to separate the neutral fields 0^ and S*^ into SU (3) octet and singlet pieces. 
The actual mass eigenstates, for Nf = 3, are mixtures which correspond more closely to the 
{uu + dd) and the 0^ (ss). 

We will work in the "heavy vector meson" limit, previously used for ordinary ChPT for 



vector mesons by Jenkins, Manohar and Wise |12]. In this approximation, the vector mesons 



are treated as static fields, whose velocity does not change when absorbing and emitting 
soft pions. The kinetic and mass terms for the vector meson fields are given by 

£km = -I ^ti{Nl {v ■ V) N^"] - I An Str[iV^] {v ■ V) Str[A^'^] (2.9) 

Cr^^s = 1^0 Str[iVt] Str[iV^] + ^ Str[iVtiV^] + Ai (Str[iVt] StT[N,M+] + h.c.) (2.10) 

+ A2Str[{iVt,iV'^}A^+]. (2.11) 

The covariant derivative is 

V^N, = d^N, + [V^,N,], (2.12) 

where 

V, = l{^d,e + ^^d,^) . (2.13) 
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There is also an axial combination, 



(2.14) 



which transforms homogeneously under the vector subgroup of the flavor symmetry. 

The propagator is somewhat different for the heavy vector mesons, but there is still a 
double pole structure in the flavor-diagonal terms. Let us neglect the SU{?>) violating mass 
terms proportional to Ai and A2. Then the propagator for the off-diagonal mesons is 



V ■ k 



(2.15) 



where = P'^ — (M^r — ^j)v^ is the "residual momentum" of the vector meson A^. For the 
flavor- diagonal fields, the propagator is 



Gt';{k) = {v>^v''-gn 



dijei ^ -Anv ■ k + fio 



V ■ k — fi 



k - /i)2 



(2.16) 



The vector mesons have a common "residual mass" fi at this order. However, we could as 
easily choose to fix /x = 0, absorbing it all into the static phase associated with the heavy 
vector meson propagator. We will do this from here on. 

There are four invariant interactions between the vector mesons and the Goldstone 
bosons: 



£int = igi Str[Art] StT[N,Ax\ v^e 



fJ,wX<T 



+ h.c. 



+ zg2 Sti[{NlN,}Ax] v^e^'^'' 
+ z(?3Str[iVt] Str[iV,] Sti[A^] v,e 
+ ig^Sii[NlK] Str[A] v^e^^^'' . 



fiuXa 



(2.17) 



Note that we could "unquench" the theory simply by removing the ghost fields and the 
flavor-singlet rj', which then would receive a large mass from the anomaly. Only the first two 
interaction terms in Ci^t would remain, corresponding to the two interactions considered in 
Ref. [|T2l. 



III. NONANALYTIC CORRECTIONS TO VECTOR MESON MASSES 

The interactions in the chiral Lagrangian effect the masses of the vector mesons in two 
ways. First, there are explicit SU{3) violating terms which correspond to local interactions in 
the theory. One cannot say much about these terms, except to use the available symmetries 
to constrain them at any given order in the chiral expansion. Second, there are "long 
distance" corrections which arise from the infrared parts of loop integrals computed using 
the Lagrangian at leading order. These corrections have a nonanalytic dependence on the 
light quark masses rrii, usually of the form ^/rr^ or rrii Inmj, and thus cannot be canceled by 
counterterms at any order in the expansion. It is useful to study these corrections, because, 
while they may not be dominant in general, they at least set a scale for SU{3) violation. As 
discussed earlier, such terms also correspond to the part of the theory where the effects of 
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FIG. 1. Feynman diagrams with no hairpin insertions. 

quenching are the least understood. Moreover, they may be dominant for sufficiently small 
quark masses. 

We will use the interaction Lagrangian >Cint to compute the nonanalytic corrections to 
the vector meson masses for Nf — 1,2,3. We will then compare our result with analogous 
corrections computed in the unquenched theory. In each case, we will truncate the expansion 
to include nonanalytic terms of order m^^, m^lnrng and m^/^, but not those of order 
rrig Inm^. Where we do not neglect them, we treat the dimensionful parameters /Iq and Mq 
as being formally of order nig. 



A. Nf = l 

The case Nf = 1 is particularly simple, since there are only flavor-singlet mesons. Since 
it is the fiavor-singlet sector which is the source of most complications in QChPT, we work 
out this case in detail as a pedagogical example. 

We begin with the propagator for the Goldstone fields. For the rj' and fj', this is given 
in matrix form by 

' '' + ^^^^(1 ])■ (3.1) 



— m^, V — 1 / (p^ — rn^'y V 1 1 

Note the relative minus sign in the non-hairpin parts of the 1]' and fj' propagators. The 
off-diagonal terms mix the r]' and the fj'. The propagator for the fermionic x fi^^d is simply 

(3.2) 



p2 — m^, 



The vector meson propagators are similar. Choosing ^ = 0, in the bosonic sector we find 
The propagators for the fermions and Kf^ are 
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FIG. 2. Quark flow diagrams which contribute to 6Ms in the quenched approximation. 



[v^'v" - g""") ^ and - {v^'v'' - g^") ^ , (3.4) 
V ■ k V ■ k 

respectively. A careful study of the ordering of creation and annihilation operators in time- 
ordered products shows that fermion loops containing an and a come with the usual 
minus sign, but those with a and a x do not. In the Feynman diagrams which we will 
consider, this effectively cancels the extra minus sign in the propagator. 

Next, we expand the interaction Lagrangian Ci^t- Since we will be computing the correc- 
tions to the real vector meson self-energy at one loop, we only keep terms with a least one 
factor of S*^. Each of the six relevant couplings has a coefficient which is a linear combination 
of the gi^s. Suppressing the common factor of v^e^'^^^ , we find the coefficients 



SlSM ■■ 2g^ + 2g2 + g3 + g4 

SlS^dxfj' : -5-3 - gi 

SlS^dxT]' : -gi - gs 

SlS^dxfj' : -gi + gs 

SlK^dxX + h-c. : gi + g2 

SlL^dxX^ + h.c. : -gi-g2. (3.5) 



There are four types of graphs which contribute to the S"^ self-energy: those with no 
hairpin insertions, those with a single hairpin insertion (on either the vector meson or the 
Goldstone boson line), and those with two hairpin insertions (one on each internal line). 
The six Feynman diagrams with only "ordinary" propagators are shown in Fig. |l|. The same 
loop integral appears in all the graphs, 

SMs oc him^i) = — -[^2^j2 ^^' (analytic in m^,) . (3.6) 

The six diagrams contribute to 6 Ms with the following coefficients of /i(m^/): 
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(a) (b) 
FIG. 3. Feynman diagrams with a hairpin insertion on the tj' propagator. 



1(a) 
1(b) 
1(c) 
E(d) 
E(e) 

0(f) 



(2^1 + 2g2 + 93 + QaY 
-{-93 -94? 
-{-91-93? 
{-9i + 93? 
-{91 + 92? 
-{-91-92? ■ 



(3.7) 



The minus signs in graphs |l|(b) and |I|(c) are from the rj' and propagators; the minus signs 
in |I](e) and [^(f ) are from fermion loops. Summing them up, we find 

{2gf + 2gl + Ag^g^ + %i^4 + %2^3 + ^929i) Hm^^,) + ... 



5Mo 



(3.8) 



Before we compute the remaining graphs, it is interesting to see which combinations of the 
giS contribute to 5 Ms- Since the role of the ghost quarks is to implement quenching by 
canceling the contributions of closed quark loops, the terms which remain after the diagrams 
have been added should be just those for which no quark loop is necessary. In Fig. |^ we 
present "quark flow" diagrams corresponding to the combinations of coupling constants 
which contribute. It is an easy exercise to check that for the combinations which do not 
appear in 5Ms-, any quark flow diagram must contain a closed quark loop. 

We pause for a moment to consider the dependence of the various terms on the number 
of colors Nc- In the large Nc limit, the coupling constants scale with Nc as follows: 



92 



1 

iv! 



1 



(3.9) 



Hence, it would be tempting to argue that the term proportional to gl in Eq. (|3.8|) is 
dominant in this limit. However, we note from Fig. |^ that the corresponding quark (and 
color) flow diagram is not planar, so it is also suppressed by l/N^.. Of course, this is what 
we would expect, since in the Nc oo limit the spectrum of the field theory is composed of 
noninteracting mesons, so Goldstone boson loops should have no effect on the vector meson 
mass. The lesson is that we must consider more than the coupling constants when extracting 
the large Nc dependence of a term in SM. 

The next class of graphs are those with a hairpin insertion on the rj' line, as shown in 
Fig. ^ In such a graph, there can be two rj' propagators, two fj' propagators, or (in two 
ways), an t]' and an fj'. So from the graph in Fig. |5(a), we find a result which is proportional 
to 



(2^1 + 2^2 + ^3 + 94? + 2(2^1 + 2^2 + ^3 + ^4) (-^3 - 94) + {-93 - 94? 



{2g^ + 2g2?, 

(3.10) 
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(a) (b) 
FIG. 4. Feynman diagrams with a hairpin insertion on the vector meson propagator. 



and from Fig. |5(b), 



-9i - 9zf + 2(-5'i - gz){-9i + Qz) + {-9i + 9zf 



-291)^ 



(3.11) 



where the minus sign in the second expression comes from the propagator. Adding these 
and including the loop integral, we obtain 



5Mc 



. . . + {4gl + ^9i92)h{mr,^) + 



(3.12) 



where 



UTTf 



(^Mqui^' - ^Aom^^ + (analytic in m^,) . (3.13) 



Similarly, we can have a single hairpin insertion on the vector meson line, as shown in 
Fig. Here the graph in Fig. §(a) is proportional to 

(2(7i + 2^72 + ^73 + 94? + 2(2(71 + 2^72 + 93 + 9a){-9i - 93) + {-9i - 93? = {9i + 2^72 + 94? , 

(3.14) 

when the and propagators are included. The graph in Fig. ^(b) is proportional to 

[-93 - 94? + 2(-5'3 - 94)i-9i + 93) + (-91 + 93?] = -{-91 - 94? , (3.15) 
where the minus sign is from the fj' propagator. Combining these terms, we find 

SMs = ... + {Agl + 4(71(72 + 4(72 (74) /a (m^') + • • • , (3.16) 

where 



127r/2 V27r 



-^iQm^i Inm^, + ANm'l^, j + (analytic in m^, 



(3.17) 



Finally, there are graphs with two hairpin insertions, one each on the r/' and 5*^ propa- 
gators. The sixteen possible graphs are illustrated compactly in Fig. |^. It is straightforward 
to enumerate them as before, and we find 



5Ms = ...+Aglh{m^, 



(3.18) 



where 
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FIG. 5. Feynman diagrams with two hairpin insertions. 



12nP 



/ioMg ( 1 + lnm^,j — —jj,QAom,^, ( 1 + 2\nm^, 



2tt 



2tx' 

3 5 



+ (analytic in rv?S) . 



(3.19) 



In Fig. P we present quark flow diagrams corresponding to the surviving diagrams with 
hairpin insertions, in which there is no closed quark loop. Again, one can convince oneself 
that all other combinations of coupling constants require quark loops, and hence do not 
contribute in the quenched theory. 



B. Af/ > 1 

If there is more than one light flavor, then the situation is complicated by the fact that 
there is more than one flavor of vector meson. However, we will see the result in this case 
can be obtained largely by considering the quark flow diagrams in Figs. || and ^ for a single 
flavor. 

The matrix of vector mesons if^ has entries, most conveniently enumerated in the 
basis qi(ly In this basis, there are three types of contributions to the mass matrix: 5Mjj, 
i ^ j, the correction to the masses of the off-diagonal meson g^g^-; 5Mjj, the correction to 
the mass of the flavor-diagonal meson q^q^] and SMu^jj, i j, the correction which mixes 
the flavor-diagonal mesons qiq^ and qjqj- Each of these corrections arises from quark flow 
diagrams with a particular structure. 

The quark flow diagrams in Figs. ^ and ^ can be divided into two classes: those in which 
the quark flavor lines run through from the incoming vector meson to the outgoing one, and 
those in which they do not. In the first class are the gig2 and g2g4 terms in Fig. |^ and the g2 
terms in Figs. |^(a) and ^(b). The others are in the second class. Because flavor is conserved 
in those diagrams in which the quark lines flow through, they contribute to 6Mij and SM^, 
but not to 6Mii-jj. In the rest of the diagrams, the quark lines in the incoming meson must 
annihilate in the graph, so they cannot contribute to the off-diagonal 6Mij] however, they 
do contribute to 6 Ma and SMu^jj. 

The only other differences in the Nf > 1 case are the somewhat more complicated 
integrals which arise when a double-pole propagator with ma ^ rrijj appears in a diagram 
with a hairpin insertion on the Goldstone boson line. The solution is written in terms of 
integrals with the following nonanalytic dependence on and y^: 
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FIG. 6. Quark flow diagrams corresponding to liairpin insertions, (a) hairpin insertion on the 
rj' line; (b) hairpin insertion on the line; (c) two hairpin insertions. 



h{x,y) 



1 



— y 
1 

1 
— 



{Ml - A,x'')h{x) - {Ml - A,y^) h{y) 
3 



Hqx lux -\- AnX 



27r 



{Ml - Aox^)h{x) - {Ml - Aoy^) h{y) 



The integrals I2 and /4 have the degenerate hmits 



h{x,x) 
h{x,x) 



1 



127r/2 
1 



Mix - -Aqx^ 



3 3 
l^oMl (1 + Inx^) - —iiqAox^ (l + 21nx^) 



27r 



27l' 

+ ^AnMIx - ^AnAqx^ 



The integrals /j will depend on the variables rriij = [BQ{mi + m^)]^/^. 

For the masses of the off-diagonal vector mesons, we find the correction 

5Mij = Agig2 h{mij) + 2g2gi [h{'mii) + h{mjj)\ 

+ 9I [h{mii, mi) + '2h{mii: rujj) + l2{mjj, rujj)] + Ag^ h{mij) . 

The flavor-conserving diagonal mass matrix is perturbed by 

5Mu = (2gl + 2gl + Agig2 + 4^i^4 + 4^2^3 + ^2^4) h{mu) 

+ (Agl + 85(15(2) l2{mii, mil) + {^qI + ^9192 + ^929^) h{mii) 
+ Aglh{mii,mii) . 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 
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Finally, the contribution to mixing between flavor-diagonal mesons is 

6Mii_jj = (^gf + 2g2g3 + ^gig^j [Ii{mii) + h{mjj)\ + 2gl h{mij) 
+ 2gig2 [hirriii, mi) + '^hirria, nijj) + hirrijj, rrijj)] 
+ (25'i5'2 + 2g2gi) [himu) + himjj)] + Agj hima, rrijj) . (3.25) 



C. Vector meson masses in ordinary chiral perturbation theory 



We would like to compare these results with the nonanalytic corrections to the vector 
meson masses computed using ordinary chiral perturbation theory. The unquenched case, 
as discussed in Ref. [0, is somewhat simpler. There are only the fields H^^ and vr, with no 
r]' field and no hairpin insertions. The interaction Lagrangian is simply 



= ig^ Tr[i7t] Tr[i7,A,] v^e^^''^'' + h.c. 



+ ig2Ti[{NlN,]Ax] v^e 



(3.26) 



Our normalization conventions are chosen to correspond as closely as possible to the 
quenched Lagrangian, and they differ slightly from those of Ref. [0]. A straightforward cal- 
culation yields the nonanalytic corrections to the vector meson mass matrix for Nf = 1, 2, 3. 
For Nf = 1, there is only the singlet S, and 



5 Ms = . 



(3.27) 



For Nf = 2, there are the isotriplet p and the singlet S. Ignoring isospin violation, we find 



6Mp = 2{gi + g2fliim^) 
5Ms = Q{gi + g2fh{m) 
6Mp.s = . 



(3.28) 



For Nf = 3, we have 

5Mk* 
SM^8 
6Ms 



2 

3' 
2 

3' 
2 

3' 
2 



3 (^1 + 1^2 ) h{m^)+gl 



3 (^1 + 1^2 ) h{mK)+gi 



3 (^1 + 1^2 ) h{mrj) + gl 



2 2 
-/i(m^) + 2Ii{mK) + -hijrir,) 

3 5 1 
-/i(m^) + -/i(mx) + -h{mri) 

2 2 
2/1 (m^) + -Ii{mK) + -himri) 



3 [gi + -92 j [3/1 (m^) + 4/1 (mi^) + /i(m^)] 
V3fi'2 {gi + ^5'2^ [3/1 (m^) - 2/1 (mx) - /i(m^) 



(3.29) 



These expressions for Nf = 3 correspond to those of Ref. [|I^ if we make the replacement 

2 \ 



9i 



1 
73 



(3.30) 
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IV. PHENOMENOLOGY 



The coefficients which appear in the chiral Lagrangians, both quenched and unquenched, 
are nonperturbative parameters which are difficult to calculate from ffist principles. How- 
ever, we need to have some estimate of their size if we are to use our expressions for phe- 
nomenology. 

The parameters Aq and Mq describing the t]' self-interaction may be studied on the lattice 
through the hairpin vertex, the anomalous scaling of m^r, and the topological susceptibility 
Xt of the pure gauge theory. The effect of the hairpin on rrij^ is parameterized by a quantity 
6, which is extracted both indirectly and directly. The values of 6 measured on the lattice 
vary widely in the approximate range < 6 < 0.3. With our normalization of the terms 
in the effective Lagrangian, the relationship between Mq and 6 is 

M^ = 8nf5, (4.1) 

from which we find < Mq < 350 MeV. 

On the other hand, Mq is also related to the topological susceptibility Xt via the Witten- 



Veneziano formula 17.18 



M^ = j^Xt. (4.2) 

One must be careful in taking this formula from the original literature, since Xt is evaluated 
in the limit Nc ^ oo. In this limit Aq = 0, that is, the wavefunction renormalization 
due to the hairpin vanishes. Hence there is no distinction between writing the left-hand 
side of Eq. as or, instead, as M^/{1 + NjAq). If we were to take this latter 

expression, then it would appear that together Eqs. (^4.1|) and ( [4.2| ) determine both Mq and 
Aq.Q However, it is inconsistent to keep the wavefunction renormalization due to Aq in one of 
these relations but not in the other. (If we do so, we are led to the contradictory conclusion 
that Aq ~ 1/-/V/, while really Aq is independent of iVj.) Instead, Eqs. ( |4.1| ) and ( |4.2| ) yield two 
independent determinations of the single variable Mq. Unfortunately, these determinations 
are not consistent with each other. Lattice calculations give Xt ~ (180 MeV)^ |TB|, from 



which we find Mq ^ 500 MeV. As a compromise, we will take Mq = 400 MeV in our 
estimates below. 

We have less information on the parameter Aq. Theoretical prejudice and the large 
Nc limit would lead one to believe that Aq is small, but perhaps the inconsistency in the 
determination of Mq is a hint that this limit does not work well here. In his recent review [115 



Sharpe fits values for a parameter = Aq/3 from a number of lattice groups, with the 
preferred value Aq ~ 0.2. This is the value we will use in the numerical estimates below, 
keeping in mind that it is probably highly uncertain. 



^In the notation of Duncan et al. [|T6|, the Lagrangian is written in terms of the parameters 
A = 1 + NfAQ and = MgNf/{l + NjAq). They take M^/{1 + NfAo) in Eq. (|]|), which 
corresponds to writing it as ttiq = ANjXt/ f'^ instead of as Aitiq = ANfXt/ f'^ ■ We agree that 
Eq. (4.1) corresponds to their expression Attiq = SNjir'^ f'^S. 
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We will also assume that / ~ 130 MeV and Bq ^ Airf are the same in quenched and 
unquenched chiral perturbation theory, and the same for the coupling constants gi and g2- 
All of the couplings Qi are expected to be of order one, with unknown sign. We will assume 
that there are no dramatic accidental cancellations between different terms, and will take 
linear combinations of the Qi^s to be of order one as well. Finally, we will neglect the vector 
meson hairpin, setting An = and yUo = 0. Our primary motivation for this is simplicity ;0 
complete expressions in which these terms are included are given above. 



A. The quenched p mass 

We now use these results to extract information about the effect of quenching in the 
determination of the p mass on the lattice. Quenched chiral perturbation theory gives 
us information about the nonanalytic dependence of rrip on the light quark masses. Let 
us consider a charged p in a theory with Nf = 2, the simplest case in which there is 
no annihilation channel. We will take the two quarks to have equal mass rriq. Then the 
expansion of rrip in terms of rUg takes the form 

mp = fi + Ci/2 + Ciml + C3/2 ml + . . . , (4.3) 

where the coefficients Ci are given by 

Ci = ~ 0.1 X (100 MeV)-^ 

-Do 



3/2 



^ Ug2igi + g4)-l0 9lAo\ ~ -0.1 x (100MeV)-^ (4.4) 



the unquenched case is C^/gm^, with 



= - ,\ '2 O.lx(lOOMeV)-^ (4.5) 



12nf 

and m^r is related to by = 2Bomq. For comparison, the leading non-analytic term in 

rrj ^ 
3/2 ''''7r5 

2ig_i±g2f 

Of course, the numerical estimates are intended only to give a sense of the order of 
magnitude of the terms. The Ci/2'^7r term is large, and its sign is predicted by the quenched 
theory. Moreover, it is the leading term for small quark mass. The other corrections are at 
the 10% level. A fit to lattice data would determine the Cj's and check whether they are 
of the expected size. In Fig. ^ we present the corrections to function of m^r, for a 

typical set of coupling constants gi of order one. We show results for two different choices 
of the couplings Mq and Aq. We also show the unquenched result. As Fig. ^ illustrates, it is 
quite generic to find extremely substantial nonanalytic corrections to rrip for ^ 400 MeV. 
Of course, the possibility of such large quenching corrections in chiral perturbation theory 
does little to enhance our confidence in the precise extraction of rrip from lattice simulations 
of quenched QCD. 



^ We note that in the real world, the size of rrip — rriuj suggests that these terms are indeed small. 
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FIG. 7. Correction 5mp to the p mass, as a function of m7r, for Nf = 2. The sohd hne is 
the result for an unquenched theory; the dashed and dotted hnes are for quenched theories with 
Mo = 0.4 GeV, = and Mq = 0.1 GeV, Aq = 0.2, respectively. We have taken all gi = 0.75. 

It is also useful to consider the case of unequal quark masses, and to study the 
quenched chiral corrections to the mass splittings. An interesting quantity to study is 
{rriK* — i^p)lijnK ~ ^V}-i which is constant in the chiral limit at leading order. The expres- 
sion simplifies if we set m,r = and study only the dependence on rriK-i in which case we 
find 

rriK* — TTip drriK* A2 



mj^ — ml dmj^ Bq 127r/^ 



mK_ 
(4.6) 



where in this expression is the one-loop corrected kaon mass. The divergence of the 
final term as m^- — > reflects the unphysical behavior of the quenched theory in the chiral 
limit. 

Finally, we may estimate the size of quenching errors by comparing the nonanalytic 
corrections to mp in the quenched and unquenched theories. Defining the ratio Rp — m^ /m^ 
to be the ratio of the masses in the two theories, we find 

Rp^l + ijl"^ [2(gf -g2 + 'ig2gi)h{m^) + ^glhim^, m^) 

- 1 - + ^Mlm^ - SAom^l , (4.7) 

for gf — 02 + 2(/25'4 = 1- For the physical meson masses m^^ = 140 MeV and = 770 MeV, 
this amounts to about a 30% error from quenching (which, of course, is dominated by the 
uncertain Mq term). On the other hand, for larger lattice masses such as ~ 500 MeV 
and mp ~ 1 GeV, the quenching error is of order 100%. Clearly, our estimate of the error 
depends on the many undetermined parameters which enter, but it does suggest that one 
ought to be cautious when using mp as a reference quantity to set the scale of the lattice 
spacing. 
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FIG. 8. The ratio rriM /rup as a function of m.,^. The vector meson parameters are as in Fig. 7, 
while the parameters for the baryon chiral Lagrangian are from Ref. [9]. 

B. The baryon- meson mass ratio 



The ratio rriN/rnp, which experimentally is 1.2, is used often to gauge the accuracy of 



lattice simulations. The chiral corrections to rrijsf have been computed by Jenkins and 
by Bernard, Kaiser and Meissner ||14|, while the corresponding quenched corrections have 
been computed by Labrenz and Sharpe P,^. Combining these with our calculation, we can 
examine the quenched chiral corrections to this ratio. We find the expansion 



= \- Ri/2 fn^ + Rim^ + Ry2 fn^ 

nip fj, 



with 



R 



1/2 



Ri 



'-{D-3Fr-Agl 



2{bD - 3bp) 



Bo 



1 



R: 



■3/2 



(D - 3F)i2D + 7) - + 4^2(^1 + 9a) 



5Ao 
3M^ 



Rl/2 



(4.8) 

(4.9) 
(4.10) 
(4.11) 



where D, F, ho-, hp, c and 7 are parameters of the baryon chiral Lagrangian |^^,13,T1|. In 
Fig. I we compare the quenched and unquenched corrections, for a typical set of coefficients 
of order one. Of course, there are far too many unknown parameters for this plot to be 
anything but illustrative of the possible size of the long distance effects. What the reader 
should note is that once again there can easily be a large difference between the two theories, 
and thus a large correction due to quenching. 
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V. DISCUSSION AND CONCLUSIONS 



Let us close by commenting briefly on our results. The most important qualitative feature 
is the dependence of nip on the square root of the quark mass rUq. This is to be contrasted 
with the unquenched theory, where the the leading dependence is linear and the leading 
nonanalytic dependence goes as m^/^. As noted earlier, while the coefficients of these terms 
depend on nonperturbative parameters, their presence is unambiguously predicted, and they 
are unchanged by the inclusion of terms of higher order in the chiral expansion. They 
dominate the extrapolation to small m^, and must be accounted for in lattice simulations. 

The derivation of our results has relied completely on chiral perturbation theory. The 
careful reader might worry that because quenched QCD is not a unitary theory, there is 
no guarantee that QChPT accurately describes its low energy limit. However, there is 
some independent evidence that QChPT is indeed the low energy theory of quenched QCD. 
The most celebrated prediction of QChPT, the emergence of "quenched" logarithms of 
the form Mq Inm^, was originally derived in another way, via strong-coupling perturbation 
theory 0,^, and the two approaches yield the same coefficient for these terms. We thus 
expect that the use of QChPT to describe quenched QCD at low energies is valid. 

We should also mention one defect of this formulation of chiral perturbation theory for 
vector mesons, namely the absence of the decay p ^ txtt. Clearly, this is an important 
contribution to the width of the p. However, its contribution to the real part of the p mass 
is rather small. Moreover, a simple one loop calculation shows that the contribution is 
nonleading, of order m^lnm-jr. We thus expect that the inclusion of this effect would not 
alter our conclusions substantially. 

Since our results depend on many unknown nonperturbative parameters, it is hard to 
draw firm quantitative lessons from them. Nonetheless, we have found in general that the 
long distance quenching corrections to vector meson masses could be quite large, at the 
level of 100%. While it is equally possible that they are small, our analysis does nothing to 
increase one's confidence that quenching is a controlled error in lattice simulations of rup and 
rriN/mp. On the contrary, the necessity of simulations of these quantities with dynamical 
fermions is as pressing as ever. 
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